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ABSTRACT 

In this Letter we examine the idea that a subset of short-period black-hole low-mass X-ray binaries could be powered by the 
mass transfer from pre-main sequence donors. As the star contracts towards the main sequence, the strong magnetic fields operate 
the magnetic braking which dissipates the orbital angular momentum, driving the binary to contact. We show that the periods and 
apparent donor spectral classes of the X-ray binaries with a pre-main sequence donor agree better with the available observations 
of black hole X-ray binaries than those of binaries with a main-sequence donor This mechanism also explains, without the need 
for additional hypotheses, the roughly primordial abundance of Li detected in donor companions of black hole X-ray candidates 
in our Galaxy. 

Subject headings: binaries; close - X-ray; binaries - stars; pre-main-sequence - stars; magnetic fields. 



1. INTRODUCTION 

Low mass X-ray binaries (LMXBs) are mass-transferring 
(MT) binaries that consist of low mass donor and a compact 
accretor, such as a black hole (BH) or a neutron star. Out of 
17 detected black-hole X-ray binaries, 1 1 have relatively short 
periods, < 15 h an d 9 of them have other properties measured 
?for reviews, see 'Casaresj 120051: iRemillard & McChntockl 
I2OO6). The effective temperatures of donor companions in 
those binaries suggest that they are low-mass stars, and their 
spectral classes indicate that they are not giants. In addi- 
tion, these short-period LMXBs are transient and in those 
of them where the lithium abundance was measured (GS 
2000H-25, A0620-00 and Nova Mus 91), observations show 
the excess above the so lar value by a factor of about 20-200 
(iMaccarone et al.ir2005L and references within). 

There are several problems with creating a complete pic- 
ture of the formation and evolution of a LMXB. The forma- 
tion scenario for such binary needs to satisfy, at the very least, 
the following obvious conditions; (i) the binary has to sur- 
vive both the supernova explosion and the pre-supernova bi- 
nary evolution; and (ii) the post-explosion binary has to be 
close enough to start the mass transfer If one requires the 
MT to start during the main sequence (MS) evolution of the 
secondary, then the latter condition limits the maximal post- 
explosion binary period to about a day (depending on the 
black hole mass); in binaries with a bigger period, the mag- 
netic braking will not provide the necessary orbital shrinkage. 

In the standard scenario for the formation of a compact 
accreting BH binary, a common envelope occurrence during 
the pre-supernova evolution is required to effect the orbital 
shrinkage. The balance between the binding energy of the en- 
velope of a massive super-giant and the released orbital en- 
ergy suggests that a low mass d onor may have diffi culties 
to eject the entire envelope (see Justham et al.l l2006b. Ap- 
pendix), unless a significant fraction of the envelope had been 
previously lost by other mechanisms; e.g., through a very ef- 
ficient stellar wind (which, on the other han d, leads to the 
binary expansion), or in an explosive merger dJustham et alj 
|2006a). Another suggested scenario, which does not have the 
energy-balance problem, is that a donor star was initially an 
intermediate mass star and the fossil magnetic fields provided 
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Strong magnetic braking that brings th e binary to the onset 
of mass transfer (lJustham et al.ll2006bh . It is crucial to un- 
derstand that the predicted formation rates for suggested sce- 
narios can vary by several orders of the magnitude, due to 
the huge uncertainties in stellar evolution (winds rates, super- 
nova kicks), the initial binary parameters (periods, mass ra- 
tios) and in the binary evolution (e.g., the common envelope 
efficiency). 

The weakness of the scenarios described above is the dif- 
ficulty to reproduce the observationally available relation be- 
tween the period and the effective temperature for the short- 
period BH X-ray bi naries. This is inevita ble if the donor of an 
intermediate mass (lJustham et al.ll20()6bh . as the theoretically 
obtained donors, being on MT sequence, are too hot compared 
to the observed ones for the same binary periods (in other 
words, for the same radii). Cold temperatures at large periods 
are also difficult to achieve for donors of smaller masses and 
occur only if the MT starts during the Hertzsprung gap, but in 
the very limited range of periods. The most difficult to explain 
are systems that have periods from 7 to 12 hours. These sys- 
tems are known to have a Lithium overabundance up to 200 
times the solar value (consistent with the primordial value). 
Note that Li overabundance has been found not in binaries 
with smallest periods within non-giant donors, but with larger 
periods, c.f. the proposition on the tidal loc king in close bi- 
naries been responsible for Li production in IMaccarone et all 
(2005) . 

We propose that these systems can be explained by consid- 
ering pre-main sequence donors. In addition to the observa- 
tional constraints on the donor temperatures and periods, we 
take into account that by the time when a primary star ex- 
plodes and creates a black hole (6-8 x 10'' yr), its low mass 
companion will not reach zero-age main sequence (ZAMS). 
In fact, at the time of 10^ yr, o nly the stars more m assive than 
2 Mq will come to their MS dP Antona & Mazzi telli 1994). 
Therefore any binary systems with a BH (and also with a neu- 
tron star) and a low-mass MS star will unavoidably evolve 
through the stage when the formed compact object is orbited 
by a pre-MS star. The radius of the star at this moment is 
about 3 times bigger than at ZAMS - the main contraction 
has already occurred. 

On Fig. [T] we show the evolutionary tracks of 1.2, 1.0 and 
0.8 Mq stars in the period — effective temperature plane, 
where the periods are such that the stars would fill their Roche 
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Fig. 1 . — Evolutionary tracks in tlie orbital period — effective temperature 
plane for donors initially of 1.2 Mq (long-dashed line), 1.0 Mq (dotted line) 
and 0.8 Mq (short-dotted Hne) from pre-MS stage; the arrows show the di- 
rection of the star's evolution through the plane. Solid line shows the position 
of ZAMS donors. All periods are such that stars would fill their Roche lobes 
if they have a 7 Mq BH companion. Dots mark the position of the donors 
at the approximate moment of BH formation, 7 X 10* yr. "Star" symbols are 
the earliest spe ctral classes of donors in observed BH X-ray binaries given in 
ICasaresI 12003) . BH X-ray binaries with Li overabundance are marked with 
circles. 

lobe with a 7 Mq companion. The evolution of stars is started 
from their Hayashi tracks. While stars evolve, they move from 
the right to the left, raising their temperature as they approach 
MS. During MS, stars perform small loops on the left from 
ZAMS line. Leaving MS, stars go from left to right and have 
much higher temperatures than at the pre-MS stage. 

It can be seen that, for the same periods, the difference 
between the hottest possible observed effective temperatures 
in Galactic BH X-ray binaries and those of the MS stars is 
very large, about 1000 K. One can argue that a star of 0.8 
Mq passes relatively close to the region of interest when it 
becomes a giant, however a star of this mass leaves the MS 
after twice the Hubble time! After the Hubble time, only 
~ 0.95 Mq star of solai" metallicity will leaves the MS. With 
a lower metallicity, Z=0.001, stars evolve faster and less mas- 
sive stars can become giants within the Hubble time. How- 
ever, it does not help, as in this case ZAMS line is shifted up- 
wards and, for the same periods, stars are hotter by about 1000 
K than in the case of Z=0.02. If a donor starts the MT during 
early MS life-time, then its evolutionary track during MT is 
going to the left and down from its MS, the track is almost 
parallel to ZAMS line on Fig. 1. This explains the observed 
BH X-ray binary with the shortest period, XTE J1180H-480 
with Prot = 4.1 h. If a star is well evolved from the MS, then 
it may follow the subgiant track and evolve to the right, but it 
will have too high temperature for the required period. When 
the star has evolved from MS, but not too far, then its track 
is similar to the behavior during the MT on MS, but it is hard 
to get far down to cold temperatures within the Hubble time. 
We must stress that, as becomes evident from the figure, the 
orbital period of a pre-MS star at the time of SN explosion has 
to be at least ^ 15h, and therefore we should eliminate from 
the consideration the periods below ~ 15h for MS stars, too. 

In what follows, we outline how the magnetic braking oper- 
ates in the pre-MS stars, show the evolution of binaries with a 
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Fig. 2. — The envelope global convective time Tc for 0.8 (open triangles) 
and 1 Mq (solid circles) stars. The dashed Hne shows the approximate time 
of the supernova explosion of the primary. Convective times are calculated 
starting from the moment when pre-main sequence stars developed a radiative 
core. Flat platos correspond to the main sequence stage. 

BH and a pre-MS star through the MT stage and analyze how 
those systems appear 

2. MAGNETIC BRAKING AT PRE-MAIN SEQUENCE STAGE 

General prediction of the dynamo theory suggests t hat ac- 
tivity should scale with the dynamo number A^d (e.g.,'ParkeiJ 
1971; Hinata 1989; Meunieret al. 1997). The simple model 
of a dynamo operating through the convective zone predicts 
that the dynamo number is related to the Rossby number 
as A^D ^ Ro~^ , where Rossby number is defined as Ro = 
Pioi/tc ijc is t he convective turno ver time and Prot is the ro- 
tation period) dNoves et al.lll984| ). The corresponding corre- 
lation between activity and the Rossby number was shown 
in a number of studies devoted to activity in late-type MS 
stars and has been used to interpret X-ray e mission of MS 
(•Montesinos et al.' '20011) and pre-MS stars (iPizzolato et al.l 
2003; Massi et al. 20 Qg), even tho ugh this correlation has a 
large scatter (Preibisch et al. ^OOS). 

On Fig.|2]we show the envelope convective turnover time, 
starting from the moment when a contracting pre-MS star de- 
velops small radiative core, and through the MS stage: is 
decreasing as the star is contracting toward the MS and re- 
mains almost the same as the star evolves through the MS. 
This figure illustrates how the correlation between Tc and 
magnetic field strength explains naturally high magnetic fields 
observed in pre-MS stars and lower magnetic fields observed 
in MS stars. As B cx v^D' ^ simple estimate shows that 
the field of a fast rotating young star with the fiot = 1 d is 
400 times stronger than that of the Sun with the Prot = 24 d: 
B^^/Bq m 400. Indeed, the magnetic field on Sun is measured 
to be 2 G, whereas t he fields detected in T Tauri stars reach 
up to 5 kG (O'Sulli vanet all 120051; ISvmington et al. 2005), 
and could be higher at earlier stages. E. g., consider TW H yd 
with measured magnetic field of 2.5 kG (lYang et alJl2005h . It 
has Teff = 4126 ± 24 K and luminosity IgL/L© = -0.62 ±0.11. 
Comparing these values with evolutionary tracks, we estimate 
TW Hyd age as ^ 2 x 10^ years. At this age t^. is several times 
less than its maximum value. 

In our study we adopt the connection between the surface 
magnetic field B^ and the Rossby number as: 
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where we take the primordial magnetic field = 10 kG; when 
dynamo mechanism starts, = 800 d and the rotation period 
is 1 d; /b is a scaling factor that allows us to vary B". This 
ratio is similar to commonly accepted relation cx il, where 
n denotes the stellar rotation rate ; used for the main sequence 
stars (e.g.. 'Mest el & Spruitlll987h . where Ro does not change 
strongly throughout the life on the main sequence. 

The magnetic braking, or the loss of the angular momentum 
from the star, is provided through the coupling of the stellar 
wind with star's magnetic field. Young stars c an lose their ma- 
terial at a very high rate, 4 x 10"^ M^lyr (e.g.. lLorenzetti et"ail 
12006). In weak T Tauri star s the mas s loss rate is as low as 
10"'"- 10"' ' MQ/yr (Andre et al.ll992h . Overall, T Tauri stars 
with s everal kG magnetic field w ould have 10"'^- 10"'' Mq/yi 
wind ( Paatz & Camenzindlll996l) . During the pre-MS stage, 
the wind strength decreases by a factor which is comparable 
to the decrease in the magnetic field strength. As the mass 
loss rate and strength of the magnetic field are correlated, 
even though not necessarily physically connected, we con- 
sider their dependence in a simple parametrized form: 

Mw = 10"'^Mo/y/w^ , (2) 

where is a parameter in the range of 0. 1 to 10. 

For the model of the magnetic braking, we adopt the stan- 
dard set of assumptions: (i) that magnetic field is dipolar and, 
accordingly, drops with distance as B(r) = B^Rr' (R is the 
radius of the star); (ii) wind corotates out to the Alfven ra- 
dius, where the magnetic pressure and ram pressure are bal- 
anced; (iii) wind-loss rate is the same at the surface and at the 
Alfven surface (for simplicity, no wind is trapped with in the 
"dead zone" like considered in lMestel & Spruij (Il987h ); (iv) 
the wind velocity at the star surface is of order of the escape 
speed. The resulting rate of change of angular momentum 
due to magnetic braking is (for the derivation see also, e.g., 
Ilustham et alJiiOOebb : 

7mb = -nB^R'^^WJ' (GM)-"' , (3) 
where M and R are the mass and the radius of the pre-MS star 

3. BINARY EVOLUTION 

We used the binary evolutionary code w ith implemented 
impli cit treatment of the mass transfer ( Ivanova & TaamI 
|2004|) . We adopted the magnetic braking described above and 
considered binaries with a 7 Mq BH with 1.0 and 0.8 Mq 
pre-MS companions. We started the evolution with several 
values of the period at the onset of MT, 9.5 h and 13 h for a 
binary with IMq companion and 13.3 h with 0.8 Mq com- 
panion. For each case, we varied the initial magnetic fields 
(/b = 1,2,5, 10 in the eq. (1)). In total, we calculated 12 evo- 
lutionary tracks. 

The results are shown on Fig [3] The end of the tracks cor- 
respond to either the stop of the mass transfer when the donor 
becomes detached; or when the donor reaches ZAMS. There 
are several reasons for the donor to become detached - it can 
contract when it approaches the MS (in the case of low mag- 
netic field), or magnetic braking stops to operate when the star 
becomes fully convective: we assume the dipole field gen- 
eration stops to be effective in this case. In the latter case. 
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Fig. 3. — Evolutionary tracks in tlie orbital period — effective temperature 
plane for donors of 0.8 Mq and 1 Mq (two case with different initial periods 
are considered in the latter case). All initial binaries are evolved with 4 values 
of magnetic pole strength (/b = 1,2,5 and 10, corresponding lines are solid, 
dotted, short-dashed and long-dashed), the resulting evolution tracks rotate 
counter-clockwise around the initial binary position (marked with large solid 
dots). Other notations are as on Fig.[T] 

the system is unlikely to reappear in the future as an X-ray 
system since the period is too wide for gravitational waves to 
shrink it, and the magnetic braking is switched off. The tracks 
shown demonstrate that it is possible to keep the binary on the 
onset of MT during the pre-MS contraction of the donor and 
that systems appear in the same region as the observed short- 
period BH X-ray binaries. 

In particular, the evolution of binaries with a low-field pre- 
MS stars is not much different from the evolutionary path they 
would have if they would be just barely at the contact and 
not losing much of their mass; they evolve to shorter periods 
and higher temperatures. As the pole-strength increases, the 
evolutionary path of the star makes counterclockwise rotation 
in the P-Titf plane, and for extremely high fields, a star is so 
much out of the thermal equilibrium that starts not only to 
cool down as it loses its mass, but also expands. 

We have estimated roughly the maximal post-explosion pe- 
riods for which a binary will start the MT during pre-MS 
stage. As expected, the higher initial B^, the wider initial bi- 
nary can be at the moment of the SN explosion and still start 
the MT during the pre-MS stage. For example, if /b = 10, 
PsN ~ 1 .8 d in a binary with 7 Mq BH and 1 Mq pre-MS star, 
and if /b = 3, then Psn is only 1.1 d. We note that those peri- 
ods are about the same or higher than those of a BH / MS star 
binary that starts MT during the MS. 

4. DISCUSSION 

We may Ust several possible criticisms of the proposed sce- 
nario: First, the interval during which such binaries can exist 
after the star formation is small, about several 10^ yr at best, 
and, accordingly, must be associated with a relatively young 
stellar population. Such a trend is not reported. On the other 
hand, if the supernova kick was as small as 50 km/s, a bi- 
nary can be several hundreds parsecs away from the place of 
it's birth. Short life-time however suggest that such systems 
can not be present among observed low mass X-ray binaries 
in elliptical galaxies and the bulg es of some spir al galaxies, 
where there is no star formation dGilfanovl 120041) . Second, 
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observed short-period LMXBs with BH donors are transient. 
Except for the systems with a low magnetic field, all other our 
systems have relatively high MT rates and should be seen as 
persistent, unless the criterion for the transiency would differ 
significantly in the case of unusually strong magnetic fields. 
And last, the life-time of the system during X-ray phase is 
shorter than it would be if a donor is at MS, about several 
10^ years versus several 10^ or even 10^ years. The short life- 
time will affect the prediction on how many binaries can be 
observed. We do not attempt to estimate the formation rates 
properly, due to the huge uncertainties that affect the rates 
for standard scenarios at first place. But we want to note that 
neither supernova kicks or the energy balance during the com- 
mon envelope would reduce the formation rates compared to 
the standard formation channel with low-mass companions, 
as they act in the same way. On the other hand, the large 
range of post-explosion binary periods may increase the for- 
mation rate. The resulting detection rates might be similar to 
those predicted for a standard channel with the same set of 
assumption for kicks and common envelope efficiency. 

In the conclusion we want to emphasize again that systems 
with a BH or a NS accretor and a pre-MS star will exist, it 
is not possible to avoid their formation. Being affected by 
strong magnetic field, those binary systems will shrink faster 
than during the MS evolution and this will alter the period dis- 
tribution of the post-supernova systems. Short period binary 
systems with a compact object and a low mass companion are 
very likely to start the MT during the pre-MS stage of the sec- 
ondary. In fact, if the post-explosion period is about a day, the 
system will definitely start the MT during the pre-MS stage 
even with a moderate magnetic field of the donor If the post- 
explosion period is larger, then the system is unlikely to start 
the MT during both pre-MS and MS stages. As one of the 
consequences, when BH-MS systems are considered, the dis- 
tribution of binary periods must take into account that pre-MS 
stars were significantly bigger at the time of the supernova ex- 
plosion. This is true as well for NS-MS systems and should 
be taken into account when, e.g., the formation of binary mil- 
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lisecond pulsars is studied. 

During the MT, the examined BH-pre-MS systems have 
the effective temperatures and periods similar to the ob- 
served BH X-ray binaries systems. As little or none nuclear 
burning occurs in donors at these ages, they will have high 
(and up to primordial) Li abundance, similar to detected in 
donors of BH X-ray binaries. It is interesting to mention 
that the don or where som e CNO-processing was seen, XTE 
J 1 1 1 8H-480 (iHaswell et al. 2002), has the period of 4. 1 h. It is 
therefore located on the left of the ZAMS line (Fig. [U and 
is explained with the evolution of a MS donor at the start 
of the MT. In addition to 3 BH X-ray binaries with Li ex- 
cess described earlier, the re are two others. One is a neu- 
ti-on star system Cen X-4 ([Martinet all 1 1994 . With the or- 
bital period of 15.1 hours, it still can be explained by our 
model. We however expect that fewer NS binaries can be 
affected by the proposed mechanism, as the formation of a 
NS occurs later than that of a BH. Another o bject is a BH 
X-ray binary V404 Cygni (iMartin et al.lll992h . With the or- 
bital periods of 6.5 days and a giant companion, this binary 
is well beyond pre-MS stage and Li excess here must have 
another origin. E.g., in low-m ass giants Li can be created via 
Cameron-Fowler mechanism (iSackmann & BoothroydlfT999l) 
with the possible enhancement in the case of tidally locked 
binaries (Denissenkov et al.l l2006h 

More detailed investigation of magnetic fields in pre-MS 
stars and their dependence on the stars age and rotation may 
help greatly for the improving of the described above sim- 
ple magnetic braking model and for the corresponding under- 
standing of LMXBs formation and evolution. We also predict 
that observations of X-ray activities of BH X-ray binaries dur- 
ing the quiescence may confirm or disproof the hypotheses on 
pre-MS nature of their donors. 
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